Dynamic symmetry breaking in a 2D electron gas with a spectral node 
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We study a 2D electron gas with a spectral node and a random gap in a vicinity of the node. 
After identifying the fundamental dynamic symmetries of this system, the spontaneous breaking of 
the latter by a Grassmann field is studied within a nonlinear sigma model approach. This allows 
us to reduce the average two-particle Green's function to a diffusion propagator with a random 
diffusion coefficient. The latter has non-degenerate saddle points and is treated by the conventional 
self-consistent Born approximation. This leads to a renormalized chemical potential and diffusion 
coefficient, where the DC conductivity increases linearly with the density of quasiparticles. 

PACS numbers: 05.60.Gg, 66.30.Fq, 05.40.-a 



The prototype of a 2D electron gas with spectral 
nodes is graphene, where two symmetric electronic 
bands, created by the underlying honeycomb lat- 
tice structure, touch each other at two different 
points in the Brillouin zone (J-Q. The surface 
states of the recently discovered topological insu- 
lators is another example for spectral nodes @. 
It is a remarkable experimental fact that the two- 
dimensional electron gas in graphene is always in a 
metallic state, regardless of its Fermi energy, pro- 
vided that the sublattice symmetry of the honey- 
comb lattice is unbroken. This is particularly sur- 
prising from the theoretical point of view, since 
the electron gas should be in a localized state, at 
least away from the node (Dirac point) |6j. At the 
Dirac point, however, the underlying Hamiltonian 
has an extra particle-hole symmetry, depending on 
the type of disorder though, which may be respon- 
sible for metallic (diffusive) behavior. The experi- 
mentally observed metallic state at and away from 
the node indicates that the diffusive behavior may 
not depend on this extra symmetry. We shall dis- 
cuss in the following that a dynamical symmetry 
exists which is responsible, regardless of the chem- 
ical potential, for a diffusive behavior. 

Many characteristic properties of a quantum sys- 
tem, in particular the spectral and large scale prop- 
erties, are determined by symmetry properties of 
the underlying Hamiltonian. However, for the dy- 
namics of a quantum system additional symme- 
tries play a role, because the dynamics is not only 
controlled by the real spectrum but also on the 
complex plane by the advanced and the retarded 
Green's function G(±ie) — (H ± ie) _1 , which 
are related by Hermitian conjugation: G(—ie) 
i • An example is the transition probability 





with 



K T>r t(ie) = (G r , r '(ie)G r /, r (-ie)), 



(1) 



(2) 



The average is here with respect to a random vari- 
able (e.g. random potential or a random gap) in 
the Hamiltonian H. Randomness is necessary to 
provide scattering that breaks translational invari- 
ance. It is convenient to combine the two Green's 
functions in the extended Green's function 



G(ie) = 



_ ((H + ie)- 









(H - ic)- 



(3) 



such that with H = diag(H, H) we have G(ie) = 
(H + iea^)^ 1 . Then there is an orthogonal or uni- 
tary symmetry which rotates the two-dimensional 
space that is spanned by the two Hamiltonians. It 
has been found long time ago that the symmetry 
breaking due to e can cause spontaneous symme- 
try breaking in the limit e — > 0. The corresponding 
massless mode leads to a diffusive behavior d, Q . 

The definition of ([3]) was also used as the start- 
ing point for Dirac fermions with random mass by 
Bocquet et al. [l(|. Employing a supersymmetric 
representation, where G(ie) is applied to a Bose 
and to a Fermi field, the gradient expansion of the 
effective field theory produces an orthosymplectic 
nonlinear sigma model in this case. Unfortunately, 
the analysis of the latter is quite involved and the 
transport properties cannot be easily extracted. 

Prior to the work by Bocquet et al., an alterna- 
tive approach was suggested by the present author 
[111 ] using explicitly the fact that the Dirac Hamil- 
tonian H — i/Jkdk+maz (cr,- are Pauli matrices and 
dj spatial differential operators) obeys the relation 



CTiiP CTi = -H 



(4) 



which constitutes class D according to Ref. [12J. 
The relation enables us to introduce the structure 



G(ie) = 



_ f(H + ieY 











(H 1 



(5) 



for the dynamic description. This choice has a 
very important advantage over ([3]), since the upper 
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FIG. 1: Typical structure of two symmetric bands 
with a node with linear (left) and quadratic (right) 
dispersion. 



and the lower block have the same determinant. 
Then the upper block can act on bosons and the 
lower block on fermions, providing us with a Bose- 
Fermi field theory and a nonlinear sigma model 
that has only a free massless Fermi (Grassmann) 
field [TH [l3[. The latter describes diffusion and 
gives directly the experimentally observed minimal 
conductivity of graphene. Moreover, it reproduces 
the phase diagram (one metal lic p hase and two in- 
sulating Hall phases) of Refs. |HH17| for a nonzero 
average mass [13| . 

The disadvantage of (JSJ over the definition ([3]) is 
its restriction to the Dirac point, since a shift by a 
chemical potential H — > H + \ioq violates the rela- 
tion (j4|). In the following we will start from ([5]) and 
extend it in such a way that a chemical potential 
can be included. This will give us a new dynamic 
structure with a continuous symmetry in Bose- 
Fermi space. The latter can be spontaneously bro- 
ken and produce a two-component massless Fermi 
(Grassmann) field. 

We consider a Hamiltonian with two bands 
whose dispersion is symmetric: ±E(k) with the 
2D wavevector k. Moreover, we assume particle- 
hole symmetry for the Hamiltonian 



OH t O 



-H 



(6) 



and include a node in the band structure (see 
Fig. [IJ. The node is important to create sponta- 
neous symmetry breaking. This has been observed 
for gapped Dirac fermions, where the symmetry- 
breaking solution vanishes when the gap is too 
large [13(. 

In analogy to the Green's function in Eq. (|5|) we 
introduce G(ie) = (H + ie)^ 1 with the extended 



Hamiltonian 



H 



Then the matrix 
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(7) 



(8) 



with scalar variables (fj, <p'j anticommutes with H: 
SH = —HS. This relation implies 



e s He s = H 



(9) 



which is a symmetry relation for the extended 
Hamiltonian with respect to U — e s . 

Now we use the matrix structure and apply it 
to a supersymmetric space that consists of four 
bosonic (upper) components and four fermionic 
(lower) components. In this representation we 
have to introduce the graded determinant detg and 
the graded trace Trg (cf. [Hj]). For constructing 
the functional integral of the transition matrix K 
in Eq. ([2j it is crucial that detg(i? + ie) = 1 and 
detg(e s ) = exp(TrgS') — 1. Moreover, the vari- 
ables ipj in (J8J) are elements of a Grassmann alge- 
bra. Therefore, the symmetry is a supersymme- 
try, connecting bosonic with fermionic degrees of 
freedom. This symmetry is broken by the e term 
though because U 2 is not a unit matrix. Eqs. ([?}, 
([8|) and (|9|) are the main results of this work. What 
remains to be discussed is the effect of the symme- 
try property on the transport for e — > 0, which will 
be studied bythe standard nonlinear sigma model 
approach (Ell, El • 

Nonlinear sigma model: The symmetry Q is 
valid for any random H or provided that H 
obeys ((6]). In order to calculate K it is necessary 
to specify the details of the randomness though. 
In the following we will consider a random gap 
because this has a strong effect due to backscat- 
tcring (unlike a random fi, where backscattering is 
suppressed). For this purpose we consider an un- 
correlated Gaussian distributed random gap with 
mean zero and variance g. 

K in Eq. ((2J can be expressed as a functional 
integral for a random term v T a a in the Hamiltonian 
H [1 31 ] (a random gap is a — 3) with the action 



S = -Trg(Q 2 ) 
9 



log detg(H + ier Q + f 3 Q) . (10) 



Q r is an 8 x 8 matrix field whose elements live 
on a superspace (the 4x4 diagonal blocks are 



3 



bosonic, the 4x4 off-diagonal blocks are fermionic) 
[Til [l3l | . We have used for the 8x8 matrix struc- 
ture, where the hat ?.. means the generalization 
of a 2 x 2 to an 8 x 8 matrix structure (e.g. 
Tj = diag(<jj,aj,<7j,<7j)), in analogy to Eq. (jTj) . 
The average Hamiltonian is Hq = (H). 

Without repeating here the derivation of the 
functional integral (cf. [l3j]) we switch directly 
to the saddle-point approximation of the integral, 
which allows us to focus on the role of the sym- 
metry in Eq. Q with the saddle point Q = irft^, 
where the shift of the chemical potential /i has been 
included in Hq as fi —¥ ft. From (|10p we obtain 

S = logdetg( J ff + i£T + t 3 Qo) , 



since the first term in (fTU|) vanishes due to 
Trg(Qo) = 0. Due to © there is a saddle-point 
manifold, when the parameters tpj are replaced 
by spatially dependent Grassmann fields spj r , that 
reads 

Qo -»■ Qr = irfi-^Ur = ivnUr • (11) 
This provides the saddle-point action 



S' = log detg ( Hq + ie + irje 



„2S 



and the transition matrix K now reads 



K , - 4— 
T 



<p jr <p' jr ,e- s V[Q] 



(12) 



(13) 



We can expand S' in powers of r\. This is also an 
expansion in powers of Go = (Ho + i(e + ti)tq) -1 , 
where the latter can be approximated by a gradient 
operator. This expansion is convergent, at least on 
large scales, when the gradient operator is small in 
comparison to r\ Up to second order in r/ it 

reads S' « S + S" with 



S" = 4i7/Trg G S 2 - 8ry 2 Trg 



GqS 



8?7 2 Trg 



(GoS*)' 



(14) 



and gives the standard form of the nonlinear sigma 
model for the off-diagonal parts of the last two 
terms, whereas the first term and the diagonal part 
of the second term contribute to the symmetry- 
breaking term that is proportional to e. Then the 
action separates into the two components as S" = 



S'2 with 



with the field $j r = ipj r (p'j r , the Laplacian d 2 = 
d\ + d 2 and the Green's function g± = [Hq + i(e + 
rf) ± ft] -1 . The parameters for e ~ read 



= -y T ^2(ff +i0 - 9-,o) = -«7 ImTr 2 (g- +fi ) , 



aj = — (— l) J a and 



(16) 



,k+q/2(^)5+,k-q/2(-«»?)] 



q=0 



(17) 

Thus our model depends only on the parameters 
g (disorder strength), r\ (scattering rate), and the 
renormalized chemical potential ft. For given g and 
IX the renormalized parameters rj and ft are deter- 
mined as a solution of the saddle-point equation. 
Interestingly, a vanishes for ft — such that the 
interaction disappears at the node (cf. [H|). Away 
from the node the interaction term in Eq. (|15[) can 
be treated in perturbation theory, where the per- 
turbations are "dimers" on the lattice which can- 
not be visited by the random walk (diffusive path). 
This restricts the random walk but since the walk 
has no phase factor (it is a classical random walk 
since (e—Dd 2 ) -1 is a real symmetric matrix), there 
is no interference to generate Anderson localiza- 
tion. 

According to ([7]) and ©, different values of j 
refer to different Fermi energies: (H + ie + /i) _1 
(for j = 1) and (H + ie - ll)- 1 (for j = 2). The 
different signs in front of a in Eq. (fl~5| reflects the 
fact that a is proportional to fi. In the following 
we ignore the index j because its value affects only 
the sign of the coupling constant . 

Random diffusion coefficient: Introducing the 
random variables u r ^ with zero mean and vari- 
ance 1, we can express the interaction term in (|15[) 
by the identity 

exp [a(5 fc $ r ) 2 ] = (exp [i\fa(d k ^r)u r ^(dk^' r )])u ■ 

which allows us to rewrite K in (|13[) . After having 
performed the ip integration we get 

K rT , = \ V ({(e-Dd 2 -i^a]u-d 2 )- 1 )) u , (19) 



9 



.7=1,2 



where the average is with respect to the new dis- 
tribution 



((...))« = (- II det(e-Dd 2 



(15) 



i.e., the determinant is included in the distribution. 
In the representation (fl~9"|) we have to average the 

r eJjU-d 2 ) -1 with 



diffusion propagator (e — Dd 2 
a random diffusion coefficient. 



This can be done in 



4 



perturbation theory by expanding the propagator 
in powers of -y/a, since a is small in the vicinity 
of the node. An alternative approach is a self- 
consistent Born approximation that formally takes 
into account a resummation of leading terms [20| • 
It is a saddle-point approximation for the integral 
over the random variables u r ^ in Eq. (|19[) . This 
approximation should be reliable, since there is no 
continuous degeneracy of saddle points as in the 
supersymmetric functional integral of Eq. (|13|) . 
The saddle-point value u implies for the diffusion 

a7u. 



coefficients Dj = D + i 
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FIG. 2: Diffusion coefficient for 2D Dirac fermions 
(in arbitray units) as a function of the renormalized 
chemical potential p for the scattering rate r\ — 0.08 
(full curve) and n = 0.1 (dashed curve). 



For a Dirac Hamiltonian (H) = i<Jkdk (valid for 
a single node in graphene or for the surface of a 
topological insulator) we get Dj = D + (— iy D' 
with 



D' ~ a\ 2 /4irD, 



where A is the momentum cut-off. The Fourier 
components of P r = P t q in Eq. (JTJ then read 



P q = KJK 



e(e + Dq 2 ) 



(20) 



(e + Dq 2 ) 2 + D' 2 q 4 
and the extension of the wavefunction ^ r r 2 P T is 



(21) 



i.e., the correction D' drops out and only the dif- 
fusion coefficient D enters the final result. D is 
plotted in Fig. [2J 

Although the dynamic conductivity is quite 
complex for Dirac particles [2l[, the DC conduc- 
tivity can be extracted from the Einstein relation 
as 



CTfcfc oc pD- 



(22) 



r?(l + r7 2 /A 2 ) ' 



with the density of states at the Fermi level p. 
Since D ~ D \p\ (cf. Fig. [2J and p ~ pn\u\ suf- 
ficiently far away from the node p = [3, |4| , the 
conductivity increases with |/2| 2 . This, on the other 
hand, is proportional to the quasiparticle density 
n. Consequently, the conductivity increases lin- 
early with n and follows the same graphs as shown 
in Fig. [21 when the chemical potential is replaced 
by the quasiparticle density. These graphs are 
in agreement with the experimental observation 

In conclusion, we have found that the transport 
properties of a 2D electron gas with a spectral node 
are controlled by a dynamic symmetry. Sponta- 
neous breaking of the symmetry generates a mass- 
less Fermi mode that leads to diffusion. For the 
special case of Dirac fermions this approach repro- 
duces the well-known V shape of the conductivity. 
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